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Abstract 

Ship  noise  measurements  are  subject  to  relatively  large  variations  that  should  be  considered 
when  these  data  are  used.  For  modern  naval  ships  operating  at  high  speed*  propeller  cavitation 
is  the  dominant  source  of  radiated  noise.  This  paper  examines  the  variability  of  the  broadband 
component  of  propeller  cavitation  noise  based  on  statistical  analysis  of  repeated  measurements 
for  several  ships  and  v  i  pdlers  and  on  a  conceptual  model  of  the  measurement  process.  Confi¬ 
dence  intervals  are  determined  for  measured  noise  spectra  and  correlations  are  sought  between 
the  variability  and  certain  important  parameters.  To  explain  these  results,  a  number  of  phe¬ 
nomena  are  discussed  that  affect  either  the  sound  transmission  properties  of  the  water,  the 
measurement  procedure,  or  the  acoustic  source  strength  of  propeller  cavitation. 

Resume 

Les  mesures  du  bruit  associe  aux  navires  peuvent  presenter  des  variations  relativement  impor- 
tantes  dont  il  faut  tenir  compte  lorsque  Ton  utilise  ces  donnees.  En  ce  qui  concerne  les  navires 
de  guerre  modernes  tres  rapides,  la  cavitation  due  aux  helices  constitue  la  principale  source  de 
bruit  rayonne.  Dans  le  present  rapport,  on  examine  la  variabilite  de  la  composante  a  large  bande 
du  bruit  associe  a  la  cavitation  due  aux  helices  on  se  fondant  sur  l’analyse  statistique  de  mesures 
repetees  avec  plusieurs  navires  et  helices  ainsi  que  sur  un  modele  theorique  de  la  methode  de 
mesure.  On  determine  des  intervalles  de  confiance  pour  des  spectres  de  bruit  mesures  et  on 
cherche  des  correlations  entre  la  variabilite  et  certains  parametres  importants.  Pour  expliquer 
les  resultats,  on  examine  un  certain  nombre  de  phenomenes  qui  influent  sur  les  caracteristiques 
de  transmission  acoustique  de  l’eau,  sur  la  methode  de  mesure  ou  encore  sur  la  force  de  la  source 
de  bruit  associe  a  la  cavitation. 
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1  Introduction 


Estimates  of  the  noise  radiated  by  naval  vessels  into  the  surrounding  water  are  required  to  eval¬ 
uate  both  the  performance  of  onboard  acoustic  sensors  and  the  probability  of  detection  by  other 
vessels.  This  information  usually  is  obtained  at  special  noise  ranges,  where  the  radiated  noise 
may  be  measured  with  the  vessel  operating  under  controlled  conditions.  Such  measurements 
often  are  required  during  the  initial  trials  of  a  new  vessel  to  enforce  the  builder’s  contractual 
obligations.  Also,  measurements  may  be  performed  at  regular  intervals  throughout  the  life  of  a 
vessel  to  evaluate  the  need  for  maintenance,  as  well  as  before  and  after  modifications  such  as  a 
change  of  propeller. 

Radiated  noise  spectra  from  sister  ships,  at  the  same  nominal  conditions,  may  be  quite  dif¬ 
ferent.  Even  the  same  ship  may  produce  noise  with  different  spectral  characteristics  at  different 
times.  This  suggests  that  ship  noise  measurements  may  include  relatively  large  random  varia¬ 
tions.  If  this  is  true,  then  any  comparison  of  ship  noise  measurements  either  with  other  data  of 
the  same  type  or  with  the  results  of  a  prediction  method  (such  as  scale-model  tests  or  computer 
models)  should  differentiate  between  variations  that  are  certain  to  be  significant  and  those  that 
may  be  due  to  unavoidable  randomness  of  the  data.  Alternatively,  if  causes  can  be  found  that 
explain  a  portion  of  this  variability,  it  may  be  possible  to  reduce  the  variability  by  controlling 
these  factors  or  by  correcting  for  their  effects. 

Figure  1  shows  a  set  of  typical  radiated  noise  spectra  for  modern  naval  surface  ships.  At  low 
ship  speeds,  groups  of  narrow  bandwidth  lines,  caused  by  various  pieces  of  rotating  machinery 
in  the  ship,  dominate  the  spectrum.  At  higher  speeds  propeller  cavitation  usually  becomes  the 
dominant  noise  source.  This  sound  has  two  distinctive  components.  Harmonics  of  the  propeller 
blade  passage  frequency  are  created  by  periodic  variations  of  the  cavitation  as  the  propeller 
turns.  Also,  broadband  sound  is  produced  with  maximum  output  at  several  hundred  Hertz.  It 
is  this  broadband  component  of  ship  noise  that  is  considered  here. 

The  goals  of  this  paper  are  to  provide  a  quantitative  understanding  of  the  variability  in 
the  broadband  component  of  full-scale  propeller  cavitation  noise  measurements,  identify  (to  the 
extent  possible)  reasons  for  the  variability,  and  provide  a  framework  for  drawing  conclusions 
from  these  data.  The  presentation  will  proceed  as  follows.  First,  Section  2  examines  the  sta¬ 
tistical  behavior  of  the  variability  based  on  analysis  of  some  currently  available  Canadian  naval 
ship  data.  Then,  Sections  3  and  4  examine  a  number  of  phenomena  which  affect  ship  noise 
measurements.  Finally,  Sections  5  and  6  make  conclusions  and  recommendations  for  good  noise 
measurement  practice. 

2  Data  Analysis 

Throughout  this  paper  the  term  ship  noise  variability  refers  to  differences  among  radiated  ship 
noise  spectra  which  cannot  be  directly  attributed  to  some  measurable  cause.  This  definition 
means  that  any  quantitative  measure  of  the  variability  will  depend  upon  how  well  the  experi¬ 
mental  conditions  are  known.  A  number  of  phenomena  are  discussed  later  in  this  paper  which 
affect  ship  noise  measurements,  but  which  are  not  typically  monitored,  and  so  are  considered  to 
contribute  to  the  noise  variability.  On  the  other  hand,  if  it  becomes  practical  to  monitor  and 
correct  for  these  phenomena,  then  the  variability  would  be  reduced. 

A  collection  of  ship  noise  spectra  was  analysed  to  estimate  the  variability.  All  of  the  data 
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Figure  1:  Typical  radiated  ship  noise  spectra  at  different  speeds 

were  measured  with  a  single  omni-directional  hydrophone  and  are  presented  as  equivalent  source 
pressure  levels  L  defined  as 

L=l0^(wi^Cr2it)PHlJ-s’)dl)-  (1) 

Here  P(tyf,Be)  is  the  acoustic  pressure  measured  at  time  t  in  a  frequency  band  of  width  Bt 
centered  at  frequency  /.  The  measurement  is  averaged  over  a  time  period  of  length  T  and  the 
pressure  is  range  corrected  from  the  actual  time  varying  range  r(t)  to  a  reference  range  ro  by 
assuming  spherical  spreading.  The  reference  range  is  chosen  to  be  one  metre.  No  corrections 
are  made  for  any  propagation  conditions  that  differ  from  spherical  spreading.  The  values  of  L 
are  expressed  in  terms  of  decibels  relative  to  a  reference  pressure  Po  of  one  micro-pascal. 

The  quantity  L  is  considered  to  be  a  random  variable  with  statistical  properties  that  are 
functions  only  of  frequency  and  propeller  revolution  rate  n.  With  these  assumptions,  the  ber 
havior  of  L  is  described  by  Pr(L  |  n,f),  the  conditional  probability  density  for  given  values  of 
frequency  and  propeller  revolution  rate  n.  The  probability  of  measuring  L  to  be  between  Li 
and  L2 ,  with  /  and  n  fixed,  is 

[ ^  Pr  (L  I  n,f)dL .  (2) 

J  L\ 

The  mean  of  L  is 

p(n,/)=  r  L?T(L\n,f)dL  (3) 

J  —  OO 

and  the  variance  of  L  is 

/) = r  pr  '  n> » dL ■  w 

J  — OO 


2 


Of  these,  the  variance  is  of  most  interest  here  because  it  is  a  measure  of  the  variability  of  the 
data. 

The  probability  density  is  estimated  below  from  experimental  data.  From  this,  the  confidence 
intervals  associated  with  single  measurements  and  with  the  mean  of  several  measurements  are 
deduced.  Then,  the  data  are  examined  for  statistically  significant  correlations  between  the 
variability  and  other  variables. 


2.1  Analysis  Procedure 


2.1.1  Classical  Estimation  Theory 


Many  random  processes  are  described  by  the  normal  probability  distribution.  If  this  is  true  for 
ship  noise  spectra  (and  later  it  will  be  shown  to  be),  classical  estimation  theory  can  be  used  in 
the  analysis  of  these  data. 

The  analysis,  as  described  in  Reference  [l],  requires  estimates  of  the  mean  and  variance  of 
the  distribution  which  can  be  obtained  from  repeated  measurements  of  L  with  fixed  n  and  /. 
If  there  are  m  repeated  measurements,  the  true  mean  p.  is  estimated  by  the  sample  mean  L, 
defined  as 


(5) 


Similarly,  the  true  variance  a2  is  estimated  by  the  sample  variance  s2, 


defined  as 


1  n 

»=i 


(6) 


Based  on  these  quantities,  the  (1  -  -y)  100  percent  confidence  interval  for  the  mean  n  is 

L  -  t^s  <  n  <  Z  +  ^Z2*  (7) 

y/m  y/m 

where  tn/ 2  is  the  value  of  the  student  t  probability  distribution,  with  m  —  1  degrees  of  freedom, 
which  encloses  an  area  of  1  —  7/2.  The  values  of  t^^fy/m  necessary  to  calculate  the  95  percent 
confidence  intervals  for  m  from  2  to  10  are  given  in  Table  1.  As  an  example  of  their  use,  the  95 
percent  confidence  interval  for  a  sample  mean  calculated  from  five  data  points  is  L±  1.24s.  This 
should  be  interpreted  as  saying  that,  if  the  set  of  five  measurements  was  repeated  100  times  to 
give  100  different  values  of  L  and  s,  then  the  intervals  L  ±  1.24s  will  contain  the  true  mean  n 
at  least  95  times  out  of  the  100. 

Part  of  the  width  of  the  above  confidence  interval  for  the  mean  is  not  associated  with 
uncertainty  in  the  data  but  with  uncertainty  in  the  estimate  of  the  variance.  The  (1  -  7)100 
percent  confidence  interval  for  the  variance  o  is 


(m  "  1)  2 

— 2 s 

Ait 


<  o' 


< 


(m-  1)  , 

— 2 - s 

*1-7/2 


(8) 


where  x*/2  and  X1-7/2  are  the  va^ues  °f  the  chi-square  probability  distribution,  with  m  —  1 
degrees  of  freedom,  which  enclose  areas  of  7/2  and  1  —  7/2,  respectively.  The  constants  needed 
to  calculate  the  95  percent  confidence  intervals  for  m  from  2  to  10  are  given  in  Table  1.  The 
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m 

to.ozs/y/m 

■smsHii 

8.98 

0.20 

1018.33 

2.48 

0.27 

39.53 

1.59 

0.32 

13.89 

5 

1.24 

0.36 

8.26 

6 

1.05 

0.39 

6.02 

7 

0.92 

0.42 

4.85 

8 

0.84 

0.44 

4.14 

9 

0.77 

0.46 

3.67 

10 

0.72 

0.47 

_ 

3.33 

Table  1:  Constants  for  95  Percent  Confidence  Interval  Calculation 


uncertainty  of  the  variance  is  very  dependent  upon  the  number  of  data  values  used  to  compute 
the  estimate  and  can  be  very  large  for  small  numbers.  The  quantity  of  data  available  for  this 
study  (typically  between  two  and  six  data  points  at  each  combination  of  frequency  and  propeller 
rotation  rate)  was  not  sufficient  to  provide  acceptable  accuracy  using  these  methods  alone. 

If  the  variance  of  the  distribution  is  known  a  priori  or  can  be  estimated  by  some  other  means, 
then  the  (1  —  7)100  percent  confidence  interval  for  the  mean  fi  is 

(9) 

where  z^/2  is  the  value  of  the  standard  normal  probability  distribution  which  encloses  an  area 
of  1  -  7/2.  The  values  of  z^j 2  are  0.67,  1.96  and  2.57  respectively  for  the  50,  95  and  99  percent 
confidence  intervals.  For  large  m,  the  student  t  distribution  asymptotically  approaches  the 
normal  distribution,  and  so  Equations  7  and  9  give  nearly  the  same  confidence  intervals.  For 
small  m,  Equation  9  gives  better  results  and  should  be  used  if  a  can  be  estimated  with  negligible, 
error.  Regression  theory  provides  a  method  to  obtain  this  improved  estimate  of  the  variance. 

2.1.2  Regression  Theory 

In  Section  2.1.1  there  was  assumed  to  be  no  correlation  between  the  statistical  properties  of  L  at 
different  values  of  the  independent  variables,  but  from  the  physics  of  broadband  cavitation  noise, 
the  mean  equivalent  source  pressure  level  should  be  a  smooth  function  of  n  and  f }  The  use  of 
regression  analysis  allows  this  property  to  be  exploited.  In  doing  so  an  assumption  must  be  made 
about  the  best  regression  model  to  use  and  the  results  depend  to  some  extent  on  this  choice. 
In  rough  terms  this  amounts  to  making  an  assumption  about  how  smooth  /i(n,  /)  is  thought  to 
be.  The  regression  analysis  methods  used  below  are  described  in  detail  in  Reference  [2]. 

Figures  1  and  2  show  the  typical  dependence  of  ship  noise  on  frequency  and  ship  speed  [3]. 
The  data  used  in  this  study  also  follow  these  trends  and  similar  plots  showed  that,  in  the 
frequency  range  above  the  spectral  peak,  the  noise  level  is  a  smooth  enough  function  of  n  and 

‘The  assumption  of  smoothness  excludes  any  consideration  of  toifal  noise  components. 
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Figure  2:  Typical  variation  of  spectrum  level  with  ship  speed  (figure  reprinted  from  Reference  [3]) 

/  that  the  following  two-variable  third-order  polynomial  regression  model  can  be  used: 

L(n,  f)  —  bo  +  f>i n  -+■  62ft2  4-  63ft3  4-  64/  4-  b$f2  4-  b$f2  4-  bjnf  -f  6gn2/  4-  b$nf2 .  (10) 

The  regression  coefficients  6o,6i,...,6g  are  determined  from  a  set  of  m  data  points  (L,-,nt,  ft) 
using  the  procedure  in  Appendix  A.  The  limitation  to  frequencies  above  the  peak  in  the  noise 
spectrum  is  not  intrinsic,  and  may  be  lifted  by  using  either  a  higher  order  regression  model  or 
this  model  in  restricted  frequency  ranges  near  the  spectral  peak  and  below  it. 

For  the  regression  analysis  the  values  L,  were  in  decibels  and  /  was  expressed  in  terms  of 
octave  band  numbers  (i.e.  proportional  to  log(/)).  This  combination  tended  to  smooth  out  the 
rapid  variations  and  allowed  the  use  of  a  low  order  regression  model. 

A  parameter  called  the  coefficient  of  multiple  determination  R2  is  a  measure  of  the  adequacy 
of  the  regression  model.  (Its  square  root  R  is  called  the  multiple  correlation  coefficient.)  It  is 
defined  as 

,  1:  (*("<,/■) -r)’ 

s^-iy 

and  has  the  value  one  for  an  exact  fit  and  zero  for  no  fit.2  Here  l!  is  the  mean  of  all  m  data 
points,  not  just  those  with  fixed  n  and  /  as  in  Section  2.1.1.  A  minimum  of  ten  data  points 
are  required  to  define  the  coefficients  in  the  regression  model.  If  only  ten  points  are  used,  the 
regression  model  will  fit  the  data  points  exactly  and  R2  =  1.  As  the  number  of  data  points  is 
increased  above  ten  the  value  of  R2  will  drop  and  will  eventually  reach  some  stable  value. 

3  A  formula  which  is  equivalent,  but  mere  efficient  for  computation,  is  given  in  Appendix  A. 
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If  the  regression  model  is  adequate  and  m  is  much  larger  than  ten,  then  L  is  an  unbiased 
estimator  of  n  with  much  smaller  confidence  intervals  than  those  found  in  Section  2.1.1.  (The 
calculation  of  the  confidence  interval  for  this  type  of  regression  model  is  quite  complex  and 
will  not  be  discussed  here.  Interested  readers  should  consult  Section  2  C  cf  Reference  [2].) 
Alternatively,  if  the  regression  model  is  not  adequate  for  some  values  of  n  or  /,  then  L  is  a 
biased  estimator.  The  adequacy  of  the  regression  model  is  determined  from  examination  of  the 
residuals  and  the  value  of  R2,  as  discussed  in  Chapter  3  of  Reference  [2]. 

Information  about  the  distribution  of  measurements  around  the  regression  value  is  found 
from  examination  of  the  residuals  which  are  defined  as 

Li- L(ntif,).  (12) 

The  variance  about  the  regression  value  is  given  by 

(is) 

t=l 


2.2  Data  and  Results 

The  analysis  procedure  described  in  Section  2.1.2  was  applied  to  actual  ship  noise  data.  The 
results  are  limited  to  frequencies  above  the  peak  in  the  cavitation  noise  spectrum  where  the 
regression  model  is  adequate.  Also,  since  broadband  propeller  cavitation  noise  wets  of  interest, 
data  were  considered  only  at  ship  speeds  where  propeller  cavitation  was  the  dominant  noise 
source  and  data  in  frequency  bands  with  strong  tonal  components  were  dropped. 

Two  sets  of  results  were  calculated.  First,  the  analysis  was  applied  to  single  data  sets  which 
are  consistent  sets  of  noise  spectra,  measured  within  a  period  of  a  few  days,  with  all  parameters 
constant  except  propeller  revolution  rate.  This  provided  information  about  the  short  term 
variability  of  ship  noise  as  would  be  encountered  during  a  typical  ship  noise  trial.  Second, 
combined  data  sets  were  formed  from  several  single  data  sets  for  the  same  ship  and  propeller, 
but  measured  at  different  times  over  a  period  of  several  years.  Analysis  of  these  provided  a 
measure  of  the  long  term  variability  which  is  most  relevant  when  comparing  data  from  several 
ship  trials. 

The  data  were  selected  from  past  ship  noise  reports  produced  by  one  of  the  Canadian  noise 
ranges  and  are  summarized  in  Table  2.  Six  combined  data  sets  were  considered,  drawn  from 
noise  measurements  of  three  ships  each  with  several  different  propeller  geometries.  For  each, 
the  propeller  cavitation  was  known  from  full-scale  trials  and  model  tests. 

2.2.1  Probability  Distributions 

The  following  procedure  was  used  to  estimate  the  probability  density  function  for  both  the 
single  and  combined  data  sets. 

1.  For  each  of  the  data  sets  with  at  least  30  points: 

(a)  apply  the  regression  model, 

(b)  calculate  the  quantity  s2,  and  then 

(c)  estimate  the  probability  distribution  of  the  residuals. 


G 


ship 

propeller 

data 

set 

m 

Single  Data  Sets 

Combined  Data  Sets 

R1 

S* 

K1 

S£ 

1 

a 

76 

0.973 

5.09 

0.955 

9.98 

b 

10 

wm 

c 

42 

0.935 

5.01 

mm 

d 

86 

0.983 

6.28 

2 

MM 

60 

0.982 

4.00 

0.937 

13.84 

MM 

51 

0.977 

6.20 

B 

1 

wm 

72 

wm 

5.46 

0.932 

9.68 

1 

36 

181 

2.11 

■9 

36 

EES] 

7.31 

2 

a 

48 

■gjkj 

8.22 

0.956 

11.36 

b 

72 

3.98 

C 

1 

a 

78 

9.48 

0.950 

9.38 

b 

12 

■I 

c 

24 

Wm 

2. 

IE 

66 

0.958 

7.88 

0.936 

11.73 

12 

n 

12 

average 

5.92 

_ 

10.99 

Table  2:  Summary  of  Data  and  Results 


2.  Form  a  composite  probability  distribution  by  combining  the  probability  distributions  found 
for  each  of  the  data  sets  and  estimate  the  overall  variance  from  the  average  of  the  separate 
values  of  s2. 

This  procedure  is  valid  only  if  both  the  probability  density  and  variance  are  independent  of  n 
and  /.  This  was  found  to  be  true,  as  discussed  in  Section  2.2.2. 

The  coefficients  of  multiple  determination  R2  and  the  estimates  of  variance  s2  are  shown  in 
Table  2.  The  smallest  value  of  R2  is  0.932  which  indicates  that  the  regression  model  is  adequate 
for  these  data. 

For  the  single  data  sets,  the  average  variance  was  5.92  (equivalent  to  a  standard  deviation  of 
2.43  dB)  and  the  composite  residual  probability  distribution  is  shown  in  Figure  3.  The  normal 
probability  distribution  with  zero  mean  and  variance  of  5.92  also  is  shown.  The  measured 
distribution  is  well  approximated  by  the  normal  distribution  especially  toward  the  tails  where 
95  percent  confidence  limits  are  located. 

The  estimated  variance  for  the  combined  data  sets  was  10.99.  The  composite  probability 
distribution  is  shown  in  Figure  4  along  with  the  normal  distribution  with  zero  mean  and  10.99 
variance.  Again  the  normal  distribution  adequately  describes  the  measured  distribution. 

These  normal  distributions  were  used  to  define  confidence  intervals  for  ship  noise  measure¬ 
ments.  The  values  for  the  50,  95,  and  99  percent  intervals  for  single  measurements  are  shown 
in  Table  3  as  given  by  Equation  9  with  m  =  1  and  a  =  s. 

The  variance  of  the  combined  data  sets  was  greater  than  that  of  the  single  data  sets.  This 
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Figure  3:  Composite  residual  probability  distribution  for  single  data  sets.  The  normal  proba¬ 
bility  distribution  with  zero  mean  and  variance  of  5.92  is  also  shown. 


Figure  4:  Composite  residual  probability  distribution  for  combined  data  sets.  The  normal 
probability  distribution  with  zero  mean  and  variance  of  10.99  is  also  shown. 
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Level 

Single  Data  Sets 

Combined  Data  Sets 

% 

(dB) 

(dB) 

50 

±1.63 

±2.22 

95 

±4.77 

±6.50 

99 

±6.25 

±8.52 

Table  3:  Confidence  intervals  for  noise  measurements  at  various  levels  of  confidence, 
suggests  that  ship  noise  variability  increases  as  the  time  between  the  measurements  increases. 
2.2.2  Correlation  of  Variability  with  Other  Parameters 

Another  objective  of  this  study  was  to  examine  statistically  significant  correlations  between  ship 
noise  variability  and  other  variables.  The  results  of  the  previous  section  provide  a  framework 
for  this. 

As  part  of  tests  for  adequacy  of  the  regression  model  that  were  described  at  the  end  of 
Section  2.1.2,  the  residuals  were  plotted  as  functions  of  noise  level,  propeller  rotation  rate,  and 
frequency.  Visual  examination  of  the  scatter  of  the  residuals  about  zero  would  show  correlations 
between  the  variability  and  these  variables,  but  no  consistent  trends  were  found  between  the 
noise  variability  and  any  of  these  variables. 

3  Noise  Transmission  and  Measurement 

A  number  of  phenomena  can  influence  ship  noise  measurements.  They  affect  either  the  charac¬ 
teristics  of 

1.  the  acoustic  source, 

2.  the  medium  through  which  the  acoustic  waves  propagate,  or 

3.  the  measurement  system. 

The  reason  for  measuring  ship  noise  is  to  determine  the  characteristics  of  the  ship  as  an  acoustic 
source.  Thus,  effects  in  the  first  of  these  categories  should  be  separated  from  those  in  the  latter 
two.  Unfortunately,  ship  noise  measurements  must  be  done  at  sea  where  it  is  seldom  possible 
to  control  or  measure  the  experimental  conditions  sufficiently  well  so  that  full  separation  is 
possible. 

In  this  section,  this  problem  is  considered  by  examining  a  conceptual  model  of  the  procedure 
commonly  used  to  measure  the  radiated  sound  from  ships.  The  purpose  of  this  analysis  is  to 
provide  realistic  estimates  of  the  errors  inherent  in  a  modern  implementation  of  the  technique. 

Figure  5  schematically  shows  the  model.  The  ship  is  assumed  to  move  in  a  straight  line  at 
constant  speed  u  past  an  omni-directional  measuring  hydrophone  located  at  depth  zm.  There 
is  assumed  to  be  a  single  compact  acoustic  source  located  on  the  moving  ship  at  depth  z,.  The 
horizontal  distance  from  the  hydrophone  to  the  closest  point  of  approach  of  the  ship  is  y  and 
the  position  along  the  ship  track  is  x.  The  pressure  at  the  hydrophone  and  ship  position  are 
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measured  as  functions  of  time.  Then  the  power  spectral  density  of  the  signal  is  computed  by 
standard  FFT  techniques  and  averaged  over  the  time  T  taken  for  the  ship  to  pass  over  the  noise 
range.  All  of  the  examples  discussed  below  have  y=100  m,  zt— 5  m,  zm=10  m,  and  u=10  m/s. 

Implicit  in  the  use  of  conventional  spectral  analysis  is  the  assumption  that  the  measured 
signal  is  statistically  stationary  and  ergodic,  which  basically  means  that  the  statistical  properties 
do  not  change  with  time.  If  this  is  true,  then  the  procedure  for  error  analysis  is  well  known.  The 
problem  here  is  that  the  signal  received  at  the  measurement  hydrophone  is  not  ergodic.  This 
is  due  primarily  to  the  motion  of  the  acoustic  source  relative  to  the  measurement  hydrophone, 
and  occurs  even  if  the  acoustic  source  itself  might  be  considered  ergodic  if  the  measurements 
could  be  done  in  a  stationary  frame  of  reference.  The  following  are  important: 

1.  The  effect  of  multiple  sound  transmission  paths  varies  with  the  position  of  the  ship. 

2.  The  measured  signal  level  changes  as  the  distance  from  the  ship  to  the  hydrophone  (range) 
changes. 

3.  Frequency  information  undergoes  a  Doppler  shift  which  is  a  function  of  the  speed  of  the 
ship  relative  to  the  measuring  hydrophone. 

The  analysis  below  deals  with  these  topics. 

3.1  The  Ideal  Measurement 

To  unambiguously  define  the  acoustic  source  characteristics,  it  is  common  to  correct  sound 
measurements  to  free-field  levels.  This  term  refers  to  the  (fictitious)  sound  measurements  that 
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would  be  made  if  the  medium  through  which  the  sound  propagates  is  both  homogeneous  and 
infinite  in  extent.  Such  a  measurement  would  include  only  sound  that  travels  in  a  straight  line 
from  the  source  to  the  receiver  and  not  any  effects  due  to  reflection,  refraction,  or  diffraction. 

If,  in  addition,  this  ideal  medium  is  lossless,  then  the  acoustic  pressure  will  always  become 
asymptotically  proportional  to  the  distance  from  the  source  after  a  certain  critical  distance  is 
exceeded.  This  type  of  pressure  variation  is  called  spherical  spreading  and  the  critical  distance 
defines  the  edge  of  the  far-field.  Near  the  source,  the  acoustic  pressure  still  may  vary  in  a 
complicated  way  even  in  an  ideal  medium. 

With  these  assumptions,  if  the  acoustic  pressure  is  known  to  be  P,(t )  at  some  reference 
distance  r o ,  then  the  acoustic  pressure  at  any  far-field  point  that  is  a  distance  r  from  the  source 

P(r,t)  = -  r/c).  (14) 

It  is  the  idealized  quantities  P,  and  its  spectral  density  S,  which  are  to  be  estimated.  Also, 
for  reasons  described  later,  the  further  condition  is  imposed  that  these  quantities  be  measured 
in  a  frame  of  reference  in  which  the  source  is  stationary. 

3.2  Measured  Pressure 

The  pressure  signal  received  at  the  measuring  hydrophone  is  a  combination  of  waves  that  travel 
by  the  direct  path  (quantities  with  subscript  d)  and  those  that  travel  by  the  surface  reflected 
path  (quantities  with  subscript  r). 

The  length  of  the  direct  path  is 

rd  =  \/ij  +  J/!  +  (z,  -  zm)2.  (15) 

Using  the  method  of  images,  the  length  of  the  surface  reflected  path  is 

rr  =  \]x7  +  V2  +  (z.  +  zmy.  (16) 

In  these  equations  and  those  to  follow,  it  is  convenient  to  introduce  the  variables 

a2  =  y2  +  (z,  -  zm)2  (17) 

and 

62  =  y2  +  (z,  +  zm)2.  (18) 

The  signal  received  at  the  hydrophone  at  time  t  is  emitted  from  the  source  at  an  earlier 
time  r  (called  the  retarded  time)  and  propagates  along  the  path  at  the  speed  of  sound  c.  In  this 
analysis  t  =  0  is  assigned  to  the  time  of  closest  approach  of  the  ship  to  the  hydrophone  with  t 
negative  as  the  ship  approaches  the  hydrophone  and  positive  after  the  ship  passes.  There  is  a 
different  retarded  time  for  each  sound  path.  The  retarded  time  for  the  direct  path  is 


Td  = 

and  for  the  surface  reflected  path 


Tr  = 


vu 


+  a2 


(19) 


t_T_Z=t_  V* 2fr  +  1)2 
c  c 


(20) 
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Here  use  has  been  made  of  the  fact  that  the  x  coordinate  of  the  ship  position  at  the  time  of 
emission  of  the  sound  is  given  by  the  ship  speed  u  times  the  appropriate  retarded  time.  Solving 
these  equations  for  the  retarded  times  in  terms  of  t  gives 


c2t  —  \/c*u2t2  +  a2(e2  —  u2) 
(c2  -  u2) 

C2t  -  \/c*U2t*  +  b2(c*  -  U *) 
(c2  -  u2) 


(21) 


The  distances  rd,  rT,  and  the  measured  ship  range  rm  may  be  expressed  in  terms  of  t  as 


Td  =  yu2rj  -f  a2 

rT  =  yj u2  r 2  -f-  62  (22) 

rm  =  \Ax2i2  +  a2. 


Here  rm  is  the  quantity  that  would  be  measured  by  an  optical  or  micro-wave  position  tracking 
system.  If  an  underwater  acoustic  range  finder  were  used,  rm  would  probably  equal  rd. 

The  measured  acoustic  pressure  is 


Pm(t)  =  Pn(t)  +  ~P.(rd)  +  Ar^P,(rr). 

rd  r. 


(23) 


This  states  that  the  measured  pressure  is  the  sum  of  waves  that  traveled  by  the  direct  path 
(second  term,  right  side),  waves  that  traveled  by  the  surface  reflected  path  (last  term,  right 
side),  and  the  quantity  P„(t)  which  includes  all  random  noise  components  such  as  the  ambient 
acoustic  noise  and  electrical  noise  in  the  instrumentation.  The  quantity  Ar  is  the  reflection 
coefficient  of  the  water  surface.  It  has  been  assumed  that  pressure  varies  with  path  length 
according  to  the  spherical  spreading  law.  Thus,  this  model  does  not  include  any  effects  of 
refraction  by  non-uniform  sound  speed  profiles. 

A  range  corrected  pressure  Pe(t)  is  then  computed  using  the  measured  range,  giving 

ft W  =  (=ft.(0  =  (=  (ft («)  +  (2ftM  +  A,?p.w) .  (24) 

ro  r0  \  rd  rT  / 

3.3  Spectra]  Analysis 


The  purpose  of  the  spectral  analysis  procedure  is  to  calculate  the  power  spectral  density  of  Pe{t). 
Currently,  the  most  common  technique  for  this  is  based  on  the  Fast  Fourier  Transform  (FFT) 
algorithm.  (See  Reference  [4]  for  a  detailed  discussion  of  FFT  analysis  techniques.) 

The  power  spectral  density  S(f)  of  the  function  Pe(t)  is  defined  as 


=  &  j 


(25) 


Here,  T  is  the  length  of  time  taken  for  the  measurement  (called  the  measurement  time).  If  Pe(t ) 
is  ergodic,  then  the  mean  may  be  converted  to  a  time  average  which  is  computed  by  dividing 
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the  fall  measurement  time  into  N  blocks  each  of  length  T'  and  starting  at  T}  .  The  equation  for 
the  estimated  power  spectral  density  then  becomes 


Se(f)  = 


2 


W  [t)Pc{t)eiwt  dt 


(26) 


The  function  W  (t)  is  a  data  window  function  which  is  introduced  to  suppress  certain  effects  that 
result  from  the  use  of  short  data  segments  [4j.  Numerous  data  windows  with  special  properties 
have  been  used  [5],  but  the  most  common  one  used  for  analysis  of  random  data  is  the  Hann 
window  defined  as 


W(t) 


=4 


1  -  cos 


(^i) 


(27) 


It  can  be  shown  that  the  operation  of  sampling  a  signal  for  a  finite  length  of  time  imposes  a 
minimum  frequency  bandwidth,  called  the  resolution  bandwidth  Be,  on  the  spectrum  calculated 
by  the  above  procedure.  This  means  that  no  feature  of  the  spectrum  can  appear  to  occupy  a 
frequency  band  less  than  Be  in  width.  The  resolution  bandwidth  is  given  approximately  by 


Bt  = 


Tr 


(28) 


This  effect  tends  to  smooth  out  sharp  features  of  the  spectrum,  and  introduces  a  bias  error  into 
the  estimated  power  spectral  density.  Reference  [4]  gives  this  as 

—  „  B\  d2S 

s‘*s+2 JaF  (29) 

To  implement  this  spectral  analysis  procedure  on  a  digital  computer,  the  function  Pc[t)  is 
first  digitally  sampled  to  produce  a  series  of  values  of  Pc(t)  at  equal  time  intervals.  This  is  called 
a  time  series.  The  integral  in  Equation  26  is  efficiently  computed  from  the  time  series  using  the 
FFT  algorithm.  The  final  power  spectral  density  is  found  by  averaging  together  a  number  of 
these  raw  spectral  estimates. 

To  avoid  a  type  of  error  called  aliasing ,  the  function  Pe(t)  must  not  contain  any  components 
at  frequencies  greater  than  half  the  sample  rate.  In  this  case  Pc(t)  is  a  measured  signal,  and 
so  it  must  be  filtered  with  an  analogue  low  pass  filter  before  it  is  digitally  sampled.  It  will  be 
assumed  that  this  has  been  done. 

Finally,  Se  is  only  a  random  estimate  of  the  power  spectral  density  which  hopefully  improves 
in  accuracy  as  more  data  is  added  to  the  average.  For  ergodic  data,  the  variance  of  Se[f)  is 


(30) 


Strictly,  this  is  only  true  if  Pc[f)  has  a  normal  probability  distribution,  but  Bendat  and  Piersol  [4] 
argue  that  it  is  a  good  estimate  no  matter  what  the  probability  distribution  as  long  as  the 
resolution  bandwidth  is  small. 


13 


3.4  Errors  in  Measured  Quantities 

Two  measured  quantities,  rm  and  Pm,  are  used  in  Equation  26  to  calculate  Se.  Each  of  these 
include  uncertainty  which  has  not  been  included  in  the  analysis  of  the  previous  section. 

A  rigorous  analysis  of  this  error  is  complicated  because  of  the  complexity  of  Equation  26, 
but  5  is  of  the  order 

'iPl  (31) 

Using  the  method  of  [6],  the  variance  of  Se  due  to  errors  in  rm  and  Pm  is 


,  ( dS  \2  f  dS  \2 

“  Ur*-;  +l aKCF-)  ■ 


where  <rrm  and  opm  are  the  variances  of  rm  and  Pm.  This  gives 


°s, 


S2. 


(32) 


(33) 


Equations  30  and  33  describe  different  components  of  the  random  error  in  Se,  but  are  not 
sufficient  to  explain  the  variance  found  in  the  statistical  analysis. 


3.5  Doppler  Shift 

The  frequency  of  the  measured  sound  differs  from  the  source  frequency  if  the  source  is  in  motion 
relative  to  the  receiver.  This  is  called  a  Doppler  shift.  The  instantaneous  measured  frequency  is 

fm  =  7~\  n  (34) 

1  -  (u/c)  CO8  0 

where  6  is  the  angle  between  the  ship  velocity  vector  and  the  direction  of  the  sound  path  at 
the  source  position  (see  Figure  6).  The  measured  frequency  is  higher  that  the  emitted  sound 
frequency  as  the  ship  approaches  the  receiver,  equal  at  the  closest  point  of  approach,  and  lower 
after  the  ship  has  passed. 

The  spectral  analysis  procedure  includes  an  average  over  the  measurement  time,  and  so  the 
Doppler  shift  is  also  averaged  over  this  time.  Figure  7  shows  the  results  for  a  1000  Hz  source 

line  with  measurement  times  of  25  and  100  seconds.  The  line  appears  spread  out.  With  very 

long  measurement  times,  9  varies  almost  from  0  to  7r,  giving  a  maximum  spread  of  14  Hz  for  a 
ship  speed  of  10  m/s.  At  shorter  times  the  spread  is  less. 


3.6  Multi-Path  Transmission 

Since  S,  is  estimated  from  ship  noise  measurements  which  must  be  made  at  sea,  it  is  important 
to  examine  how  the  ocean  differs  from  a  true  free-field  environment  and  determine  how  this 
affects  the  results. 

Sound  is  refracted  in  the  ocean  due  to  variations  in  water  temperature,  pressure,  and  salin¬ 
ity.  In  some  circumstances  this  may  affect  ship  noise  measurements,  but  here  it  is  assumed 
that  the  distance  from  acoustic  source  to  the  measurement  hydrophone  is  small  enough  that 
refractive  effects  can  be  ignored.  If  the  distance  is  less  than  one  kilometer,  this  should  be  a  good 
assumption. 
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-HYDROPHONE 


Figure  6:  Sound  path  geometry 


Figure  7:  Effect  of  Doppler  shift  on  a  1000  Hz  source  line  for  measurement  times  of  25  and  100 
seconds 
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The  reflection  of  sound  does  affect  ship  noise  measurements  significantly.  The  air-water 
surface  is  a  good  reflector  of  sound,  and  so  surface  reflections  are  usually  important  and  are 
treated  in  the  following  discussion.  Reflections  from  the  sea  bottom  may  also  be  important, 
but  are  more  difficult  to  treat.  Also,  if  the  measurements  are  carried  out  in  deep  water  or  over 
areas  of  sea-bed  with  high  acoustic  absorption,  reflections  from  the  sea  bottom  may  be  negligible. 
Thus,  bottom  reflections  are  not  treated  in  detail,  but  their  effects  should  be  qualitatively  similar 
to  those  of  the  surface  reflections. 

Figure  6  shows  the  geometry  of  the  sound  paths.  Sound  that  travels  by  the  surface  reflected 
path  takes  slightly  longer  to  reach  the  hydrophone  than  the  direct  path  sound.  The  resulting 
phase  difference  causes  an  interference  effect  sometimes  called  the  Lloyd’s  mirror  effect. 

As  shown  in  Appendix  B,  if  the  source  and  receiver  are  both  stationary,  the  power  spectral 
density  of  the  measured  sound  ST(f)  is  related  to  S,(f)  by 


Sr(f)  = 


1  +  2 Ar  —  COS  (Jjfjd  -  Tr)  +  A2/-^ 

rr  rr 


(35) 


Here  the  acoustic  source  is  moving,  but  the  spectral  analysis  procedure  treats  the  time  series  in 
short  segments  over  which  the  difference  in  time  delays  is  nearly  constant.  Thus,  the  measured 
power  spectral  density  Sm  is  approximately  the  average  of  Sr  over  the  total  measurement  time, 
and  so 


Sm(f)  _  1 

s,(f)  T 


ta  r2 

1  +  2 Ar  —  cos  u(rd  -  rr)  +  A]  -f 

Tr  rr 


dt. 


(36) 


The  quantity  Ar  is  the  reflection  coefficient.  For  a  flat  water  surface  it  has  a  value  very  close  to 
minus  one  for  all  angles  of  incidence,  but  for  rough  surfaces  Ar  becomes  a  function  of  angle  of 
incidence  and  frequency.  As  shown  in  Reference  (7],  a  good  approximation  is 


Here,  a,  is  the  root  mean  squared  surface  roughness  and  f)  is  the  angle  of  incidence  (see  Figure  6). 
The  value  of  a,  is  approximately  the  significant  wave  height  divided  by  four.  The  significant 
wave  height  can  be  found  from  sea  state  tables  [8],  Sea  state  zero  would  correspond  to  a  at  of 
zero,  for  sea  state  one  the  maximum  a,  is  approximately  0.025  m,  and  for  sea  state  two  it  is 
0.125  m. 

Equation  36  gives  the  error  in  the  measured  power  spectral  density  that  is  induced  by  the 
effect  of  the  surface  reflection.  The  integral  was  evaluated  numerically  and  the  results  for  a 
flat  sea  surface  (reflection  coefficient  of  minus  one)  and  no  attenuation  are  shown  in  Figure  8 
as  a  function  of  frequency  for  several  measurement  times.  At  low  frequencies  the  error  is  a 
strong  function  of  frequency  and  the  length  of  time  taken  for  the  measurement.  While  at  high 
frequencies  the  error  gradually  steadies  down  to  nearly  3  dB. 

This  error  is  also  a  function  of  sea  state.  Figures  9  and  10  compare  the  error  due  to  surface 
reflections  for  a  flat  sea  to  that  with  a  a,  of  0.125  m  (top  of  sea  state  2).  Figure  9  was  calculated 
with  a  measurement  time  of  25  seconds  and  Figure  10  with  100  seconds.  In  both  cases  the 
rough  sea  surface  reduces  the  high  frequency  reflection  coefficient  so  that  error  due  to  surface 
reflections  is  virtually  eliminated  at  50  kHz.  The  results  with  the  longer  measurement  time  show 
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Figure  8:  Lloyd’s  mirror  correction  vs.  frequency  for  measurement  times  of  25,  50,  75,  and  100 
seconds 

more  error  at  high  frequencies  because  more  data  is  included  in  the  average  from  times  when 
the  source  is  far  away  and  f)  is  large  (resulting  in  a  large  reflection  coefficient). 

These  results  were  calculated  without  attenuation,  but  its  inclusion  would  not  greatly  affect 
this  component  of  error  because  the  attenuation  is  nearly  equal  for  the  two  paths.  Also,  Equa¬ 
tion  36  cannot  describe  a  Doppler  shift.  This  would  be  important  at  high  frequencies  where  it 
would  smooth  out  the  ripples  in  the  spectrum. 

3.7  Attenuation 

Another  important  factor  in  the  measurement  of  source  strength  is  the  attenuation  that  occurs 
as  the  sound  propagates  along  its  path.  The  pressure  varies  as 

P  =  P,  exp(-aer).  (38) 

Here  ci!e  is  the  attenuation  coefficient.  It  is  more  common  to  express  this  in  terms  of  the 
logarithmic  attenuation  coefficient  a  which  is  found  from 

^20  log  =  8.686ae  =  o,  (39) 

where  a  is  expressed  in  units  of  dB  loss  per  unit  distance. 

Near  the  ocean  surface  there  are  two  important  components  of  this  attenuation:  the  atten¬ 
uation  of  pure  sea  water,  and  that  due  to  suspended  bubbles. 
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Figure  9:  Lloyd’s  mirror  correction  vs.  frequency  for  rms  sea  roughnesses  of  0  and  0.125  m, 
measurement  time  is  25  seconds 
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Figure  10:  Lloyd’s  mirror  correction  vs.  frequency  for  rms  sea  roughnesses  of  0  and  0.125  m, 
measurement  time  is  100  seconds 
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Figure  11:  Attenuation  coefficient  vs.  frequency  for  various  water  temperatures  (salinity  of  35 
parts  per  thousand) 

3.7.1  Pure  Sea  Water 

The  attenuation  of  pure  sea  water  is  a  function  of  frequency,  temperature  and  salinity.  The 
variation  with  frequency  at  various  temperatures  is  shown  in  Figure  11,  as  calculated  by  the 
method  in  Reference  [7]  for  a  salinity  of  35  parts  per  thousand  by  weight.  The  attenuation 
of  pure  sea  water  is  very  small  below  5  kHz,  increases  roughly  as  /2  to  nearly  100  kHz,  and 
decreases  with  increased  temperature. 

The  results  of  including  this  attenuation  in  the  error  calculation  of  Equation  36  are  shown 
in  Figures  12  and  13.  Figure  12  shows  the  results  for  calm  water,  temperatures  of  5  and  20  C, 
and  a  25  second  measurement  time.  The  attenuation  is  greatest  at  high  frequency  and  the 
difference  in  temperature  causes  approximately  0.8  dB  difference  in  error  at  50  kHz.  The  longer 
measurement  time,  shown  in  Figure  13,  causes  the  effect  to  be  more  pronounced  because  the 
average  then  includes  sound  that  travels  over  longer  paths. 

3.7.2  Suspended  Bubbles 

The  sound  attenuation  just  described  applies  deep  in  the  ocean.  Near  the  surface  the  water  may 
contain  small  quantities  air  bubbles  which  greatly  increase  the  attenuation  over  that  of  pure  sea 
water. 

A  number  of  measurements  of  the  number  and  size  of  bubbles  in  sea  water  are  reviewed 
in  Reference  [9].  Near  the  water  surface  the  action  of  wind  and  wa'es  cause  small  air  bubbles 
to  be  mixed  into  the  water  column.  The  large  bubbles  rapidly  rise  to  the  surface  leaving  a 
distribution  of  smaller  bubbles  which  rise  less  quickly.  In  deep  water  far  from  land  this  is  the 
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Figure  12:  Lloyd’s  mirror  correction  vs.  frequency  for  water  temperatures  of  5  and  20  C, 
measurement  time  is  25  seconds 


Figure  13:  Lloyd’s  mirror  correction  \s.  frequency  for  water  temperatures  of  5  and  20  C, 
measurement  time  is  100  seconds 
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dominant  method  of  bubble  formation  and  the  number  of  bubbles  per  unit  volume  (bubble 
density)  varies  exponentially  with  water  depth  and  is  a  function  of  wind  speed.  In  shallow 
water  near  land  the  bubble  density  does  not  diminish  as  quickly  with  depth  nor  is  it  as  strong 
a  function  of  wind  speed,  suggesting  other  sources  of  bubbles  such  as:  wave  activity  against  the 
coast,  biological  activity,  or  bubbles  carried  into  the  sea  from  rivers. 

The  motion  of  a  surface  ship  also  adds  bubbles  to  the  water.  The  wake  behind  the  ship  is 
an  area  of  turbulent  water  containing  many  bubbles  which  may  persist  for  hours  after  passage 
of  the  ship.  Also,  the  rolling  motion  of  the  ship  may  trap  air  under  portions  of  the  hull  such 
as  the  bilge  keels.  Thus,  a  surface  ship  which  repeatedly  passes  through  the  same  area  of  water 
may  increase  the  bubble  density  in  that  area. 

The  attenuation  caused  by  these  bubbles  is  very  dependent  on  both  the  sound  frequency  and 
the  bubble  sizes.  Laboratory  measurements  with  water  containing  bubbles  all  of  nearly  equal 
size  show  that  the  attenuation  near  the  resonance  frequency  of  the  bubbles  is  over  100  dB  per 
meter  for  air  contents  as  low  as  0.05  percent  by  volume,  but  drops  off  quickly  at  frequencies  away 
from  resonance  [10,11].  Sea  water,  on  the  other  hand,  will  contain  bubbles  of  many  sizes  which 
resonate  at  different  frequencies,  and  so  will  have  a  relatively  broad  peak  in  the  attenuation 
curve. 

The  attenuation  that  is  actually  measured  in  any  particular  circumstance  cannot  be  predicted 
unless  the  distribution  of  bubble  sizes  along  the  path  the  sound  takes  from  the  source  to  the 
measurement  point  is  known.  Novarini  and  Bruno  [12]  attempted  such  a  calculation.  They 
derived  an  empirical  formula  for  the  bubble  size  distribution  for  the  deep  ocean  as  a  function  of 
water  depth  and  wind  speed,  which  they  used  to  study  refraction  and  attenuation  of  sound  by 
the  surface  bubble  layer.  They  showed  the  bubbles  could  form  a  sound  channel  at  the  surface 
and  estimated  the  attenuation  per  limiting  ray  cycle  (approximately  4  km)  of  10  kHz  sound 
trapped  in  this  layer  as  nearly  6  dB  for  15  knot  winds  and  between  35  to  40  dB  for  25  knot 
winds.  The  attenuation  of  sound  that  leaves  the  bubble  layer  was  found  to  be  strong  function 
of  the  source  and  receiver  position  in  the  layer. 

These  results  show  such  a  strong  dependence  on  the  sound  path  and  bubble  size  distribu¬ 
tion  that  practical  calculations  are  very  difficult.  Still,  there  is  good  reason  to  believe  that 
surface  bubbles  significantly  affect  ship  noise  measurements  if  either  the  acoustic  source  or  the 
hydrophone  are  within  the  bubble  layer.  The  amount  of  attenuation  may  vary  significantly  over 
a  short  period  of  time  (several  hours)  in  response  to  changes  in  the  bubble  distribution  brought 
about  by  changing  wind  conditions  or  the  build  up  of  bubbles  introduced  by  the  ship  wake. 

3.8  Combined  Errors 

So  far  each  type  of  error  has  been  considered  separately,  but  in  reality  they  are  all  combined. 
Figure  14  compares  the  combined  errors  due  to  surface  reflections  and  pure  sea  water  attenuation 
for  two  different  (but  reasonable)  situations.  One  case  is  for  a  flat  sea,  20  C  water  temperature, 
and  a  25  second  measurement  time.  The  other  is  for  sea  state  2,  5  C  water  temperature,  and  a 
100  second  measurement  time.  The  ship  noise  spectra  measured  under  these  conditions  would 
differ  by  as  much  as  7  dB  at  some  frequencies.  The  effect  of  suspended  bubbles  was  not  included 
here,  but  would  lead  to  even  greater  variability. 
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Figure  14:  Error  calculated  for  two  specific  cases 

3.9  Choice  of  Measurement  Time 

Previous  discussions  have  shown  that  the  measurement  time  is  an  important  parameter,  but  did 
not  indicate  how  to  choose  it. 

The  most  basic  consideration  is  that  the  measurement  time  be  long  enough  so  that  the 
uncertainty  of  the  spectral  density  (given  by  Equations  9  and  30)  meets  the  desired  limit. 
Also,  if  the  ship  is  rolling,  the  noise  measurement  should  extend  over  several  roll  periods  (each 
on  the  order  of  10  seconds).  If  this  amount  of  data  cannot  be  obtained  in  a  single  pass  by  the 
hydrophone,  then  spectra  from  additional  passes  should  be  averaged  until  the  total  measurement 
time  meets  the  requirement. 

All  of  the  examples  of  errors  due  to  surface  reflections  and  attenuation  were  calculated  with 
a  ship  speed  of  10  m/s.  At  other  speeds,  the  same  error  curves  apply  as  long  as  the  geometry  of 
the  sound  paths  is  the  same.  This  requires  that  the  noise  be  measured  over  a  constant  portion 
of  the  ship’s  track  past  the  hydrophone.  (In  other  words  the  noise  should  be  averaged  over  a 
constant  range  of  x  coordinates.)  This  requires  the  measurement  time  to  vary  inversely  with 
ship  speed. 

4  Propeller  Cavitation 

So  far  effects  related  to  the  transmission  and  measurement  of  sound  have  been  examined,  but 
not  variation  of  the  actual  acoustic  source  level.  This  section  deals  with  factors  which  affect 
propeller  cavitation  as  a  noise  source. 

The  cavitation  can  change  in  two  ways.  First,  the  ship  speed  at  which  it  first  occurs  (called 
the  inception  speed)  may  change.  Second,  at  speeds  above  inception,  new  types  of  cavitation 
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may  form  or  the  relative  quantity  of  each  of  the  types  may  vary. 

These  changes  in  the  cavitation  can  be  produced  in  several  ways. 

1.  If  the  propeller  is  controllable  pitch,  the  pitch  must  be  accurately  set,  which  may  be 
difficult  to  do  accurately. 

2.  Any  damage  to  the  propeller  or  fouling  (perhaps  by  marine  growth,  rope  or  fishing  nets) 
would  greatly  increase  the  cavitation. 

3.  A  change  in  the  depth  of  the  propeller  (for  example  caused  by  a  change  in  ship  draught) 
alters  the  cavitation  inception  speeds. 

4.  A  change  in  the  propeller  inflow  velocity  distribution  will  vary  both  the  types  of  cavitation 
and  the  inception  speeds.  This  could  be  caused  by  modifications  to  the  hull  especially  the 
addition  of  appendages,  by  cross  flow  introduced  by  ship  motions  or  water  current,  or  by 
rudder  action  required  to  maintain  course. 

5.  An  increase  in  the  ship  resistance  will  cause  more  cavitation  due  to  the  increased  loading 
of  the  propeller. 

If  the  geometry  of  the  ship  and  propeller  are  constant,  then  the  last  two  of  these  factors  are 
most  important.  Of  these,  it  is  possible  to  quantify  the  effect  of  ship  .resistance  on  cavitation. 

4.1  Ship  Resistance 

As  discussed  above,  the  ship  resistance  is  an  important  factor  in  determining  the  amount  of 
propeller  cavitation  and  thereby  the  noise  level.  To  illustrate  this,  an  estimate  of  the  change  in 
noise  due  to  a  10  percent  increase  in  ship  resistance  is  presented. 

Suppose  that  back  sheet  cavitation  is  the  dominant  form  of  caviitation.  The  discussion  in 
Reference  [13]  (Section  16.7  of  Chapter  VII)  suggests  that  a  10  percent  increase  in  ship  resistance 
would  increase  the  inception  cavitation  number  by  20  percent.  The  shape  of  the  noise  level  versus 
cavitation  number  curve  is  required  to  estimate  the  resulting  increase  in  noise.  Reference  [14] 
includes  such  a  curve  for  hydrofoil  surface  cavitation,  which  shows  that  an  increase  in  inception 
cavitation  number  of  this  size  could  cause  an  increase  in  the  sound  levels  of  over  10  dB. 

This  calculation  should  be  considered  as  only  a  rough  estimate.  Better  results  could  be 
achieved  by  considering  model  test  results  for  particular  ships,  but  the  conclusion  would  still 
be  that  changes  in  ship  resistance  of  the  order  considered  here  can  significantly  affect  the  noise 
output  from  propeller  cavitation.  Several  factors  that  affect  the  ship  resistance  will  now  be 
discussed. 

4.1.1  Shallow  Water 

In  very  deep  water,  ship  performance  is  independent  of  the  actual  depth.  This  is  not  true  in 
shallow  water.  The  nearby  ocean  bottom  alters  both  the  surface  wave  pattern  produced  by  the 
ship  and  the  velocity  of  the  water  beneath  the  hull.  This  causes  the  draught,  trim,  and  ship 
resistance  to  change  from  the  values  in  deep  water. 
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Figure  15:  Percentage  increase  in  ship  resistance  in  shallow  water  (figure  reprinted  from  Refer¬ 
ence  [13]) 

Detailed  methods  of  estimating  the  change  in  ship  resistance  in  shallow  water  are  given  in 
Reference  [13]  and  Figure  15  illustrates  a  specific  example.  A  critical  parameter  is  the  dimen¬ 
sionless  quantity 

Fn  =  TfiK’  (40) 

where  u  is  the  ship  speed,  g  is  the  acceleration  due  to  gravity,  and  h  is  the  water  depth.  The 
effect  of  water  depth  is  negligible  if  this  quantity  is  less  than  0.4.  The  ship  resistance  reaches  a 
maximum  at  values  of  Fn  between  0.8  and  1. 

As  an  example,  suppose  that  a  ship  with  5  metre  draught  is  in  30  metres  of  water.  Curve  B 
of  Figure  15  suggests  that  the  ship  resistance  exceeds  the  deep  water  value  by  10  percent  when 
Fn  is  greater  than  0.7,  which  corresponds  to  a  ship  speed  of  12  m/s  (23  knots). 

4.1.2  Other  Effects 

During  the  normal  operation  of  a  ship  the  draught  and  trim  may  change  due  to  consumption 
of  fuel  or  addition  of  equipment  and  cargo.  These  change  the  wetted  surface  area  and  possibly 
the  surface  wave  pattern,  which  will  affect  the  ship  resistance. 

The  resistance  also  depends  upon  the  condition  of  the  hull  surface.  A  ship  is  usually 
smoothest  and  has  the  lowest  resistance  immediately  after  the  hull  is  painted  and  increases 
after  that.  For  conventional  paints  the  increase  in  resistance  may  be  from  20  to  30  percent  in 
the  first  year  of  operation  [13].  For  newer  self-polishing  anti-fouling  paints  this  is  not  such  a 
problem. 
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4.2  Environmental  Factors 

It  is  best  to  measure  ship  noise  in  flat  seas  with  no  wind  or  water  current.  Unfortunately,  these 
ideal  conditions  are  seldom  available.  The  following  discussion  considers  how  poorer  conditions 
affect  ship  noise. 

Water  current  does  not  affect  ship  performance  other  than  to  change  the  velocity  vector.  A 
current  directed  along  or  opposite  the  ship  heading  alters  the  speed  over  the  ground,  but  not 
the  radiated  noise.  On  the  other  hand,  a  cross  current  may  cause  greater  propeller  cavitation  if 
significant  rudder  must  be  applied  to  keep  the  ship  on  the  required  course. 

Compared  to  this,  wind  usually  has  a  much  larger  effect.  The  force  applied  to  the  ship  is 
proportional  to  the  squared  magnitude  of  the  apparent  wind  (wind  measured  relative  to  the  ship 
velocity  vector)  and  varies  with  direction.  For  a  ship  moving  directly  upwind  or  downwind,  this 
force  changes  mainly  the  ship  resistance.  Reference  [13]  suggests  that  an  apparent  head  wind 
of  15  m/s  will  increase  resistance  from  5  to  20  percent  (depending  on  the  superstructure)  over 
that  with  zero  apparent  wind.  Wind  that  is  not  directly  on  the  bow  or  stern  will  cause  the  ship 
to  list  and  yaw,  as  well  as  require  significant  rudder  to  maintain  course. 

Finally,  a  rough  sea  influences  all  aspects  of  ship  performance  and  greatly  increases  propeller 
cavitation.  The  ship  resistance  increases,  as  do  all  of  the  components  of  ship  motion  (such  as 
roll,  pitch,  etc.).  These  effects  cause  the  propeller  loading  to  vary  as  the  ship  moves  in  the  sea, 
which  causes  the  amount  of  cavitation  to  vary.  Also,  air  can  be  trapped  under  parts  of  the  hull, 
such  as  bilge  keels,  and  be  swept  into  the  propeller  causing  great  changes  in  the  noise  produced 
by  the  cavitation. 

5  Conclusions 

This  paper  examined  the  variability  of  far-field  noise  from  ship  propeller  cavitation.  An  analysis 
of  some  naval  ship  data  provided  statistical  information  about  the  variability,  but  revealed  little 
about  its  causes.  Further  discussion  of  phenomena  that  affect  ship  noise  measurements  showed 
that  the  statistical  results  were  reasonable  and  identified  a  number  of  areas  where  ship  noise 
measurement  practice  could  be  improved. 

The  normal  probability  density  function  adequately  described  the  variability  of  the  noise 
spectra  considered  here.  The  variability  did  not  depend  upon  noise  level,  frequency,  or  propeller 
revolution  rate,  but  the  time  between  measurements  was  important.  Data  sets  measured  over  a 
period  of  years  had  greater  variability  than  sets  containing  data  taken  over  a  few  days. 

The  confidence  intervals  in  Table  3  are  based  on  these  results.  They  are  valid  when  broadband 
cavitation  noise  is  the  dominant  noise  source,  for  frequencies  above  the  spectral  peak.  These 
confidence  intervals  provide  a  good  indication  of  the  variability,  but  are  based  only  on  data  for 
three  ships  measured  at  one  noise  range.  More  general  results  would  require  data  from  more 
ships  and  noise  ranges. 

Part  of  the  variability  is  due  to  random  error  in  the  spectral  analysis  procedure,  and  uncer¬ 
tainty  in  the  measurements  of  sound  pressure  and  distance  from  the  ship  to  the  hydrophone. 
These  are  an  inherent  part  of  any  noise  measurement  and  are  well  understood.  Unfortunately, 
they  are  insufficient  to  explain  the  variability  of  ship  noise  data. 

The  sound  transmission  properties  of  the  ocean  cause  part  of  the  additional  variability.  Sea 
roughness  affects  the  amount  of  high  frequency  sound  reflected  from  the  ocean  surface.  Water 
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temperature  and  salinity,  as  well  as  the  number  of  suspended  gas  bubbles,  affect  the  attenuation 
of  sound.  These  phenomena  can  cause  ship  noise  measurements  to  vary  by  more  than  5  dB, 
sufficient  to  explain  the  statistical  results.  This  could  happen  over  a  period  of  several  hours. 
Thus,  sound  transmission  measurements  would  be  required  several  times  a  day  to  monitor  these 
effects. 

The  geometry  of  the  sound  paths,  which  is  determined  by  the  positions  of  the  ship  and 
hydrophone,  also  affects  the  transmission  properties.  This  source  of  variability  can  be  minimized 
by  always  measuring  the  noise  over  a  constant  portion  of  the  ship’s  track  past  the  hydrophone. 

Finally,  the  noise  produced  by  the  cavitation  may  vary.  This  may  result  from  changes  in  ship 
resistance  due  to  operation  of  the  ship  in  shallow  water,  as  well  as  changes  in  hull  roughness, 
draught,  and  trim.  Also,  deviation  from  the  ideal  environmental  conditions  of  flat  seas  with  no 
wind  or  water  current  will  cause  greater  noise. 

6  Recommendations  for  Ship  Noise  Measurement 

The  following  recommendations  will  help  control  the  variability  of  ship  noise  measurements  and 
should  be  implemented: 

1.  Noise  ranges  should  calculate  confidence  intervals  for  their  measurements,  possibly  using 
the  procedure  described  in  Section  2.1.2.  The  distinction  between  short  term  and  long 
term  variability  is  important  and  separate  confidence  intervals  should  be  determined  for 
each. 

2.  The  sound  transmission  properties  of  the  water  should  be  measured  using  broadband 
acoustic  sources  located  at  several  positions  along  the  ship  track  and  at  the  same  depth 
as  the  ship’s  dominant  noise  source.  This  should  be  done  several  times  a  day  especially  if 
the  weather  changes. 

3.  The  test  program  should  be  done  in  random  order  to  prevent  changes  in  sound  transmission 
properties  from  introducing  systematic  errors. 

4.  The  noise  should  be  measured  over  a  constant  portion  of  the  ship’s  track  past  the  hy¬ 
drophone. 

5.  In  all  but  flat  seas,  the  noise  measurement  should  include  data  from  several  ship  roll 
periods  (each  around  10  seconds). 

6.  Both  shaft  power  and  revolution  rate  should  be  measured  to  indicate  changes  in  ship 
resistance. 

7.  Environmental  conditions  such  as  sea  state,  wind,  and  water  temperature  should  be 
recorded. 

8.  The  maximum  ship  speed  should  be  less  than  0.4 \/gH.  At  speeds  above  this  limit,  ship 
resistance  depends  strongly  on  water  depth,  causing  the  cavitation  noise  to  be  different 
from  that  in  deep  water. 
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9.  Water  currents  which  cause  the  ship  to  apply  significant  rudder  to  maintain  course  should 
be  avoided.  On  the  other  hand,  currents  along  or  opposite  the  ship  heading  do  not  affect 
noise.  Currents  in  any  direction  sweep  away  air  bubbles  introduced  into  the  sound  range 
area  by  the  ship’s  motion. 
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Appendix  A 

Given  a  set  of  m  data  points  of  the  form  (ii,y«,2t),  a  solution  is  sought  for  the  regression 
coefficients  in  the  regression  model 

z  =  b0  +  bxx  +  62i2  +  63is  +  64  y  +  bsy2  +  bey3  +  67IJ/  +  b&x2y  +  bgxy2 .  (A.l) 

The  regression  coefficients  are  found,  as  shown  in  Reference  [2],  from  the  solution  of  the  matrix 
equation 

AB  =  C.  (A.2) 

For  the  stated  regression  model 
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The  solution  procedure  consisted  of  first  constructing  the  matrices  A  and  C  from  the  m  data 
points,  and  then  determining  the  matrix  of  regression  coefficients,  B,  using  a  numerical  subrou¬ 
tine  for  solution  of  simultaneous  linear  equations. 
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The  coefficient  of  multiple  determination  is 


,  mBTC-G»» 
where  Br  is  the  transpose  of  matrix  B. 


(A.6) 
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Appendix  B 

In  this  appendix,  the  power  spectral  density  of  a  signal  that  contains  a  reflection  is  calculated. 
The  signal  is  assumed  to  be  of  the  form 

f{t)  =  fo{t)  +  Afo{t-t0),  (B.l) 

where  /o(t)  is  the  unreflected  signal,  A  is  the  strength  of  the  reflected  signal,  and  to  is  the 
amount  the  reflected  signal  is  delayed. 

The  Fourier  transform  of  this  is 

F(u)  =  I*  fo(t)eiutdt  +  A  [°°  /0(t  -  t0)eiutdt.  (B.2) 

J  —  oo  J  — oo 

Changing  variables  to  t'  =  t  -  to  in  the  second  integral  gives 

F(u)  =  fo{t)eiuidt  +  Aeiut°  f°°  fo(f)eiut' dt' .  (B.3) 

J — oo  J  —oo 

Both  of  these  integrals  are  the  Fourier  transform  of  /o(t),  and  so 

F(v)  =  (l  +  Ae‘"‘0)  F0(v).  (B.4) 

The  power  spectral  density  is  proportional  to  |F(w)|2  which  is 
|F(w)|2  =  F(u,)F*  (w) 

=  (l  +  Ae^0)  (l  +  Ae-’“‘0)  |FoH|2  (B.5) 

=  [l  +  A  («**•  +  e~iut°)  +  A2]  |F0(u;)|2. 

Now,  Euler’s  Formula  can  be  used  to  replace  the  exponentials  with  2  cos(wto)  giving 

|F(w)|2  =  ^1  +  2Acos(wto)  +  A2]  |F0(w)|2.  (B.6) 
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